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Abstract 

Using an isospin- and momentum-dependent transport model we study near-threshold pion pro- 
duction in heavy-ion collisions induced by radioactive beams at the planned Rare Isotope Accel- 
erator (RIA). We revisit the question of probing the high density behavior of nuclear symmetry 
energy E sym (p) using the 7r~ /tt + ratio. It is found that both the total and differential ir~ /tt + 
ratios remain sensitive to the E sym {p) when the momentum-dependence of both the isoscalar and 
isovector potentials are consistently taken into account. Moreover, the multiplicity and spectrum 
of 7r~ mesons are found more sensitive to the E sym {p) than those of ir + mesons. Finally, effects of 
the Coulomb potential on the pion spectra and 7r~ /ir + ratio are also discussed. 
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I. INTRODUCTION 



The isospin asymmetry of pions produced in nuclear reactions was shown to be useful for 
extracting interesting information about the structure of radioactive nuclei and neutron skins 
of heavy rtahie nuclei within the Gianher .nodeiflflflflB. Within hanronic transport 
models the isospin asymmetry of pions was found useful also for studying the equation of 
state (EOS) of isospin asymmetric nuclear matter^, Q, Isj]. In particular, the 7t~/tt + ratio 
was proposed as a sensitive probe of the high density behavior of nuclear symmetry energy 
-Esj/m(p)|Z|- The latter is an important part of the nucleon specific energy E(p, 8) 

E(p, 5) = E{p, 5 = 0) + E sym (p)S 2 + 0(S 4 ) (1) 

in asymmetric matter of isospin asymmetry 5 = (p n — Pp)/(p p + p n )- The density 
dependence of symmetry energy is very important for many interesting astrophysical 



problems 
ion reactions 



ice oi symmetry ( 

300,0,13, 



ry interesting a; 
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the structure of radioactive nuclei \14. Il5l Il6l 1171 1 and heavy- 
4|, LLi| • Unfortunately, the density dependence of symmetry energy, especially 
at supranormal densities, is still very poorly known. Predictions based on various many- 
body theories diverge widely at both low and high densities. In fact, even the sign of 
the symmetry energy above 3po is still very uncertain[12|. Since nuclear reactions induced 
by high energy radioactive beams can produce transiently dense neutron-rich matter, the 
fast fragmentation beams from RIA and the new accelerator facility at GSI provide the 
first opportunity in terrestrial laboratories to explore experimentally the EOS of dense 
neutron-rich matter 4], Q, 18, 0|- Crucial to the extraction of critical information about 
the E sym (p) is to compare experimental data with transport model calculations. Based 
on isospin-dependent transport model calculations, several experimental observables have 



been identified as promising probes of the E, 
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V), such as, the neutron/proton ratio 



isoscaling in nuclear mult if r agment at ion \2 It |22j, 23], the neutron-proton differential flow 



Q, 



the neutron-proton correlation function |28( and the isobaric yield ratios of 
light clusters[29]. The first experimental constraint on the density dependence of symme- 
try energy at subnormal densities was recently obtained by analyzing the isospin diffusion 



data from heavy-ion reactions 
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3l|. While it is much more challenging to constrain the 



symmetry energy at supranormal densities. In anticipation of coming experiments with high 
energy radioactive beams at RIA and GSI, more theoretical studies identifying experimental 
observables sensitive to the symmetry energy at supranormal densities are needed. 



The single nucleon potential is an important input to the transport models. It includes an 
isovector part (symmetry potential) and an isoscalar part, and both of them are momentum 
dependent due to the non-locality of strong interactions and the Pauli exchange effects 
in many-fermion systems. However, in all transport models for heavy-ion collisions the 
momentum-dependence of the symmetry potential was seldom taken into account until very 
recently. It was found that the momentum dependence of the symmetry potential affects 
manv experimental observables that were known to be sensitive to the symmetry energy (31L 
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33J. In this work we study near-threshold pion production at RIA using a momentum and 



isospin dependent transport model 32] . We revisit the question of whether the high density 
behavior of nuclear symmetry energy E sym (p) can be probed using the it~/tt + ratio. We 
found that both the total and differential n~ /tt + ratios remain sensitive to the E sym (p) when 
the momentum-dependence of both the isoscalar and isovector potentials are consistently 
taken into account. Moreover, the multiplicity and spectrum of tt~ mesons are found more 
sensitive to the E sym (p) than those of 7r + mesons. 



II. A BRIEF SUMMARY OF THE ISOSPIN AND MOMENTUM DEPENDENT 
TRANSPORT MODEL IBUU04 

In this work, we use the isospin and momentum dependent transport model for heavy-ion 
collisions induced by neutron- rich nucleiS. In the latest version of this model IBUU04 we 



use a single nucleon potential 



U(p,5,p,T,x) = A u (x) — + Ai(x) — 

Po Po 



+ ) CT (1 - x5 2 ) - 8tx - P<J 1 5p T , 
Po ° + 1 Po 



po J "l + (p-^) 2 /A 2 

Po J l + (p-^) 2 /A 2 
In the above r = 1/2 (—1/2) for neutrons (protons) and r ^ r'; a = 4/3; f T (r,p) is 
the phase space distribution function at coordinate r and momentum p. The parameters 
A u (x), Ai(x), B, C t ,t> CV,t' and A were obtained by fitting the momentum-dependence of 
the U(p, S,p, t, x) predicted by the Gogny Hartree-Fock and/or the Brueckner-Hartree-Fock 
calculations, saturation properties of symmetric nuclear matter and the symmetry energy 



of 30 MeV at normal nuclear matter density po — 0.16//m 3 3|. The compressibility of 
symmetric nuclear matter Kq is set to be 211 MeV. 

The last two terms contain the momentum-dependence of the single particle potential. 
The momentum dependence of the symmetry potential stems from the different interaction 
strength parameters C T y and C TiT for a nucleon of isospin r interacting, respectively, with 
unlike and like nucleons in the background fields. More specifically, we use C un iike = —103.4 
MeV and Cuke = —11.7 MeV. One characteristic of the momentum dependence of the 
symmetry potential is the different effective masses for neutrons and protons in isospin 
asymmetric nuclear matter. With the above potential, we found that the neutron effective 
mass is higher than the proton effective mass and the splitting between them increases with 
both the density and isospin asymmetry of the mediumj^]. Moreover, with the potential of 
eqUJboth the isoscalar and isovector potentials at po are in agreement with the corresponding 
nucleon optical potentials extracted from nucleon- nucleus scattering data 3^1 • 

The parameters A u (x) and Ai(x) depend on the x parameter according to 

IB 

AJx) = -95.98 -x , (3) 

a + 1 

2B 

A t {x) = -120.57 + x -. (4) 

IT + 1 

The parameter x can be adjusted to mimic predictions on the density dependence of the 
symmetry energy E sym (p) by microscopic and/or phenomenological many-body theories. In 
this work we choose a range of the x parameter from 1 to -2. With this choice the density 
dependence of the symmetry energy samples a wide range of theoretical predictions as shown 
in Fig. 1. 

Other details of the IBUU04 model can be found in ref. 32j|. The details of modeling 

nn 

pion production can be found in our earlier publications, e.g., ref . |7|, |35[ . Pion production 
was also studied previously by many other people within transport models. However, none 
of the previous studies has taken into account the momentum dependence of the symmetry 
potential. Only very recently the latter was recognized as an important issue in connection 
with probing the high density behavior of symmetry energy using the 7r~ /tt + ratio. 
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III. RESULTS AND DISCUSSIONS 



We study the reaction of 132 Sn + 124 Sn at a beam energy of 400 MeV/nucleon and an 
impact parameter of 1 fm as an example of high energy central reactions at RIA. This partic- 
ular reaction will be carried out using the fast fragmentation beam line at RIA. The beam 
energy selected here is about the highest one to be available at RIA. A Time Projection 
Chamber (TPC) has been proposed to study the EOS of isospin asymmetric nuclear matter 
by detecting charged particles including pions[36(. In the following we discuss several fea- 
tures of near-threshold pion production at RIA. The emphasis will be on examining their 
sensitivities to the density dependence of symmetry energy. 

A. Formation of high density isospin asymmetric nuclear matter at RIA 

What are the maximum baryon density and isospin asymmetry that can be achieved in 
central heavy-ion collisions at RIA? This is an interesting question relevant to the study of 
the EOS of asymmetric nuclear matter. To answer this question we show in Fig. 2 the central 
baryon density (upper window) and the average (n/p) p > Po ratio (lower window) of all regions 
with baryon densities higher than p . It is seen that the maximum baryon density is about 2 
times normal nuclear matter density. Moreover, the compression is rather insensitive to the 
symmetry energy because the latter is relatively small compared to the EOS of symmetric 
matter around this density. The high density phase lasts for about 15 fm/c from 5 to 20 
fm/c for this reaction. It is interesting to see that the isospin asymmetry of the high density 
region is quite sensitive to the symmetry energy. The soft (e.g., x = 1) symmetry energy 
leads to a significantly higher value of {n/p) p > po than the stiff one (e.g., x = —2). This is 
consistent with the well-known isospin fractionation phenomenon. Because of the E sym (p)5 2 
term in the EOS of asymmetric nuclear matter, it is energetically more favorable to have 
a higher isospin asymmetry 5 in the high density region with a softer symmetry energy 
functional E sym (p). In the supranormal density region, as shown in Fig. 1, the symmetry 
energy changes from being soft to stiff when the parameter x varies from 1 to -2. Thus 
the value of {n/p) p > Po becomes lower as the parameter x changes from 1 to -2. It is worth 
mentioning that the initial value of the quantity {n/p) p > Po is about 1.4 which is less than 
the average n/p ratio of 1.56 of the reaction system. This is because of the neutron-skins of 
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the colliding nuclei, especially that of the projectile 132 Sn. In the neutron-rich nuclei, the 
n/p ratio on the low-density surface is much higher than that in their interior. It is clearly 
seen that the dense region can become either neutron-richer or neutron-poorer with respect 
to the initial state depending on the symmetry energy functional E sym (p) used. 

B. Dynamics of pion productions in heavy-ion collisions induced by radioactive 
beams 

To understand the dynamics of pion production and its dependence on the symmetry 
energy, we show in Fig. 3 the multiplicity of tt + , tt~ and A (1232) as a function of time. The 
multiplicity of A(1232) resonances shown in the figure includes all four charge states while 
in the model we do treat and follow separately different charge states of the A (1232) and 
iV*(1440) resonances. At a beam energy of 400 MeV/nucleon which is just about 100 MeV 
above the pion production threshold in nucleon-nucleon scatterings, almost all pions are 
produced through the decay of A(1232) resonances. The contribution due to N* resonances 
is negligible. By comparing Fig. 2 and Fig. 3 one can notice that most of the A resonances 
are produced in the high density region. Pions from the decay of these resonances thus 
carry useful information about the high density phase. It is interesting to see that the ir~ 
multiplicity depends more sensitively on the symmetry energy. This feature is seen more 
clearly in Fig. 4 where the n~ and tt + multiplicities at the freeze-out are shown as a function 
of the x parameter. While the 7r + multiplicity remains about the same the ir~ multiplicity 
increases by about 20% by varying the x parameter from -2 to 1. The multiplicity of tt~ is 
about 2 to 3 times that of 7r + . This is because the 7r~ mesons are mostly produced from 
neutron-neutron collisions. The ir~ mesons are thus more sensitive to the isospin asymmetry 
of the reaction system and the symmetry energy. In fact, assuming pions are all produced 
through A(1232) resonances in the first chance nucleon-nucleon scatterings and neglecting 
the influence of subsequent, piou scattering* aud reabsorptions, the ratio is expected 

to scale with the N/Z ratio of the participant region according to|37| 

7T-/7T+ = (5iV 2 + NZ)/(5Z 2 + NZ) « (N/Z) 2 . (5) 

In reactions induced by neutron-rich nuclei one thus expects to see more 7r~ than 7r + mesons. 
For the reaction considered here, the average (N/Z) 2 of the reaction system is about 2.4. 
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The observed 7r~ /n + ratio is somewhat different from this value depending on the symmetry 
energy used. It indicates the effects of the different N/Z ratios of the participant regions 
with the different x parameters due to the isospin fractionation. Moreover, it also indicates 
the importance of pion reabsorption and rescatterings which are both taken into account 
properly in our transport model investigations. 

Our finding that ir~ mesons are more sensitive to the symmetry energy has interesting 
implications for the experimental studies. With the TPC, for instance, the ir~ mesons can 
be easily identified because of their negative charges. While the 7r + mesons will be bent to 
the same direction in the magnetic field as protons and are thus difficulty to be separated 
out cleanly. To investigate further the possibility of exploring the symmetry energy using 
7r~ mesons alone, we show in Fig. 5 the 7r~ and 7r + kinetic energy spectra. It is seen that the 
7r~ spectra with different x parameters can indeed be clearly separated. It thus raises the 
interesting possibility of probing the density dependence of symmetry energy by studying 
the 7T~ spectrum. 



C. The 7r /ir + ratio probe of the high density behavior of symmetry energy 



We now turn to the 7t~/tt + ratio as a probe of the high density behavior of symme- 
try energy. The advantage of using the 7r~ /h + ratio over the 7r~ spectrum itself is that 
the ratio can reduce largely the systematic errors involved in the experiments. Moreover, 
within the statistical model for pion production 0, the 7r ~/ 7r+ ratio is proportional 
to exp [(/!„ — fi p )/T], where T is the temperature, \i n and /i p are the chemical potentials 
of neutrons and protons, respectively. At modestly high temperatures (T > 4 MeV), the 
difference in the neutron and proton chemical potentials can be written as 4(| 



Pn Pp 



jjn _ jjv 

^ asy w asy 



Coulomb 



+ T 



ln^ + Y ^B n 



(6) 



where Ucouiomb is the Coulomb potential for protons, At is the thermal wavelength of a 



40|. The 



2U asy 5, where the 



nucleon and B' m s are the inversion coefficients of the Fermi distribution function 
difference in neutron and proton symmetry potentials U™ sy — U^ sy 
function U asy is the strength of the symmetry potential, depends on the x parameter, the 
density p and the nucleon momentum. It can be estimated readily from the single nucleon 
potential of eq|21 It is also seen that the kinetic part of the difference \i n — fi p relates directly 



to the isospin asymmetry p n j p v or p n — p p . Thus from the statistical point of view the 
7r~/7r + ratio is theoretically a good probe of the symmetry energy. However, to verify this 
expectation and study more realistically the ratio one has to rely on the transport 

models. 

Shown in Fig. 6 is the quantity /^ + )uke a s a function of time. Taking into account 
the dynamics of resonance production and decays we define 

, *- + A- + jA° ( 7) 
v 1 ! 7T+ + A++ + |A+ y ! 

This ratio naturally becomes the final tc~ /ir + ratio after all resonances have decayed. First, it 
is seen that the {n^ /7i + )uke ratio reaches a very high value in the early stage of the reaction. 
This is due to the abundant neutron-neutron scatterings when the two neutron-skins start 
overlapping at the beginning of the reaction. Secondly, the (vr _ /7r + ) Kfce ratio saturates after 
about 25 fm/c indicating that a chemical freeze-out stage has been reached. Finally, the 
sensitivity to the symmetry energy is clearly shown in the final tt~ /tt + ratio. One can notice 
that the sensitivity of the 7r~ /h + ratio to the x parameter is quantitatively about the same 
as the 7T~ multiplicity shown in Fig. 4. Shown in Fig. 7 and Fig. 8 are the differential tt~ /n + 
ratios versus the kinetic energy and transverse momentum, respectively. In the low energy 
{Ekin < 120 MeV) or transverse momentum (p t < 200 MeV/c) region the n~ /n + ratio is 
clearly separable with the x parameter varying from 1 to -2. It indicates that it is sufficient 
to measure accurately the low energy (transverse momentum) pions, instead of the whole 
spectrum, in order to constrain the density dependence of the symmetry energy. 



D. Coulomb effects on the it — /tt+ ratio 
First of all, it is necessary to mention that the Coulomb effects on tt~/tt + ratio is wel 
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known, in particular from experiments at the BEVALAC and SIS/GSI, see, e.g., ref. 
and references therein. A number of models have been used in analysing the n~ /n + ratio in 
heavy-ion reactions from low to ultra-relativistic energies. However, in previous studies the 
effects of the symmetry energy was generally neglected. From eq. E]it is seen that the 7r~/7r + 
ratio depends on both the Coulomb and symmetry potentials. Here we are interested in 
understanding the relative effects of the symmetry and Coulomb potentials. As an example, 
we compare in Fig. 9 the multiplicities of tc~ and tt + mesons as a function of transverse 



S 



momentum calculated with and without the Coulomb potential with the parameter x — 0. 
It is seen that the Coulomb potential is to shift n~ to lower (higher) p t as one expects. 
The 7r~ //T+ ratio with and without the Coulomb potential is shown in Fig. 10. It is seen that 
the ir~ /tt + ratio is quite flat without the Coulomb potential. The Coulomb effect increases 
the 7i~/'K + ratio at p t = for about 30% while suppresses it at higher p t . What is the 
relative effect of the symmetry potential with respect to that of the Coulomb potential? 
The answer to this question can be obtained from inspecting Fig. 7 and Fig. 8. With 
different x parameters, the reaction dynamics is about the same as indicated by the almost 
identical evolution of central density shown in Fig. 2. Thus the Coulomb effect should also 
be about the same in calculations with the different x parameters. From the variation of the 
7r~/7r + ratio by varying the x parameter as shown in Figs. 7 and 8, we can conclude that the 
Coulomb effect is stronger than the symmetry potential. Of course, we should stress that 
the symmetry potential has an indirect effect on the n~ /n + ratio through its interactions on 
nucleons, while the Coulomb potential acts directly on charged pions. Since the Coulomb 
effect is well-known, our results discussed above indicate that it is very promising to extract 
useful information about the symmetry potential from studying the 7i^/n + ratio. 



IV. SUMMARY 



In summary, using an isospin- and momentum- dependent transport model we have stud- 
ied near-threshold pion production in heavy-ion collisions induced by radioactive beams at 
the planned Rare Isotope Accelerator (RIA). We studied properties of the high density mat- 
ter formed in high energy central reactions at RIA. It is found that the isospin asymmetry 
of the high density hadronic matter is very sensitive to the symmetry energy used. With 
the soft symmetry energy the n/p ratio of the high density region can be significantly higher 
than that of the reaction system due to the isospin fractionation. We also examined the 
multiplicities and spectra of both n~ and n + mesons. We found that the n~ mesons carry 
more sensitive information about the symmetry energy than the 7r + mesons. Both the ki- 
netic energy and transverse momentum spectra of 7r~ mesons can be useful for studying 
the symmetry energy at high densities. Moreover, we revisited the question of probing the 
high density behavior of nuclear symmetry energy E sym (p) using the 7i~ /tt + ratio. It was 
found that both the total and differential n~ /ti + ratios remain sensitive to the E sym (p) when 
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the momentum-dependence of both the isoscalar and isovector potentials are consistently 
taken into account. Furthermore, we found that the effects of the Coulomb potential on 
the tt~/tt + ratio are important. In fact, the Coulomb effect is stronger than that of the 
symmetry potential. Nevertheless, taking into account the Coulomb effect and the momen- 
tum dependence of the symmetry potential the 7t~/tt + ratio, especially in the low energy 
or transverse momentum region, remains very sensitive to the variation of the symmetry 
energy. Our study thus confirms that the isospin asymmetry of pions from heavy-ion reac- 
tions induced by radioactive beams is a very promising tool for studying the EOS of dense 
neutron-rich matter, especially the high density behavior of symmetry energy. Our findings 
here are expected to be useful for designing the TPC and planning experiments at RIA. 
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FIG. 1: (Color on line) Nuclear symmetry energy as a function of density with different x parameter. 
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FIG. 2: (Color on line) Central baryon density (upper window) and isospin asymmetry (lower win- 
dow) of high density region for the reaction of 132 Sri + 124 Sn at a beam energy of 400 MeV/nucleon 
and an impact parameter of 1 fm. 
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FIG. 3: (Color on line) Evolution of the pion and A(1232) multiplicity in the reaction of 132 Sn+ 124 
Sn at a beam energy of 400 MeV/nucleon and an impact parameter of 1 fm. 
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FIG. 4: (Color on line) The average multiplicity of ir + and ir~ as a function of the x parameter 
for the reaction of 132 Sn + 124 Sn at a beam energy of 400 MeV/nucleon and an impact parameter 
of 1 fm. 
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FIG. 5: (Color on line) The kinetic energy spectra of it and ir + in the reaction of 132 Sn + 124 Sn 
at a beam energy of 400 MeV/nucleon and an impact parameter of 1 fm. 
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FIG. 6: (Color on line) Evolution of the 7r~ /tt + ratio in the reaction of 132 Sn + 124 Sn at a beam 
energy of 400 MeV/nucleon and an impact parameter of 1 fm. 
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FIG. 7: (Color on line) The ir /tt + ratio as a function of pion kinetic energy in the reaction of 
132 Sn + 124 Sn at a beam energy of 400 MeV/nucleon and an impact parameter of 1 fm. 
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FIG. 8: (Color on line) The tt /tt + ratio as a function of transverse momentum in the reaction of 
132 Sn + 124 Sn at a beam energy of 400 MeV/nucleon and an impact parameter of 1 fm. 



20 



t 1 r 




p t (MeV/c) cm. 



FIG. 9: (Color on line) ir~ and tt + transverse momentum spectra calculated with and without the 
Coulomb potential for the reaction of 132 Sn + 124 Sn at a beam energy of 400 MeV/nucleon and 
an impact parameter of 1 fm. 
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FIG. 10: (Color on line) 7r~ /tt + ratios as a function of transverse momentum calculated with 
and without the Coulomb potential in the reaction of 132 Sn + 124 Sn at a beam energy of 400 
MeV/nucleon and an impact parameter of 1 fm. 
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